Composite pulses in NMR quantum computation 
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I describe the use of techniques based on composite rotations to combat systematic errors in 
quantum logic gates. Although developed and described within the context of Nuclear Magnetic 
Resonance (NMR) quantum computing these sequences should be applicable to other implementa- 
tions of quantum computation. 



INTRODUCTION 

Quantum computers [H, 0| are explicitly quantum me- 
chanical systems that use quantum phenomena to per- 
form computational tasks more efficiently than any clas- 
sical computer. Unsurprisingly quantum computation 
has generated enormous interest, reflecting not just its 
potential technological importance, but also the intellec- 
tual importance of the challenge provided to previous 
formulations of computational complexity theory. This 
interest is tempered by the apparent difficulty in building 
large scale devices capable of implementing useful com- 
putations, but it has proved fairly simple to build small 
demonstration devices, and nuclear magnetic resonance 
(NMR) has played a leading role in this. 

The first ideas on how to build quantum computers 
with NMR 0, Si, 0,13 were swiftly followed by the first 
implementations of quantum algorithms [1, i, [l^, [HI . It 
must be remembered that the great difficulty in prepar- 
ing NMR systems in pure spin states has given rise to 
grave concerns about the direct relevance of NMR tech- 
niques to attempts to build large scale devices 12, l^, 1^ . 
and has even led to questioning of whether many NMR 
quantum computations can be considered true quantum 
computations at all [l5[. Despite this, NMR remains 



an interesting technology for exploring simple quantum 
phenomena and developing techniques which may find 
applications in other technologies. A number of reviews, 

e.g. [B, [13, HE US], have described developments in 

particular areas of NMR quantum computing. In this 
paper I will describe the use of composite pulses, a tech- 
nique developed in conventional NMR 2]J , to the design 
of quantum logic gates which are resistant to systematic 
errors in their implementation 2^ 13, 3, 25| . 



SPIN-HALF NUCLEI IN LIQUID SAMPLES 

There are many possible physical implementations of a 
qubit, but a particularly natural implementation is pro- 
vided by a spin- half atomic nucleus, such as ^H. A two- 
qubit quantum computer can be built from two atomic 
nuclei, and so on. It is necessary that the two nuclei are 



distinct, so that the the two qubits can be separately ad- 
dressed, and there must be some sort of spin-spin inter- 
action, so that two-qubit logic gates can be constructed. 
This is easily achieved by using two inequivalcnt nuclei 
in a molecule. 

The spin Hamiltonian [2€ 
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[29| is in principle 
quite complicated, but in liquid state (or solution state) 
samples is greatly simplified by rapid molecular tum- 
bling. This largely removes intcrmolecular interactions 
such as dipolar coupling, and cancels out the anisotropic 
parts of intramolecular interactions, such as the chemical 
shift and scalar coupling, reducing them to their isotropic 
forms. The cancellation of intcrmolecular interactions re- 
sults in an ensemble of identical independent molecules, 
which can for most practicalpurposes be treated as a sin- 
gle molecule in a (usually [30|) highly mixed spin state. 
(This cancellation is not, of course, perfect, and inter- 
molecular and anisotropic interactions remain an impor- 
tant source of spin relaxation.) 

An important practical distinction can be made be- 
tween systems where all the spins are of different nuclear 
species (a fully heteronuclear spin system) and those with 
two or more nuclei of the same type (a homonuclear spin 
system). As there are a limited number of different spin- 
half nuclei, among which only six (^H, ^^C, ^^N, ^^F, ^^Si 
and ^^P) have been used in quantum computing exper- 
iments so far, it is clear than only very small quantum 
computers can be fully heteronuclear, but the relative 
ease of working with such systems makes them popular 
for implementing simple tasks. 

Systems involving quadrupolar nuclei, with spins 
greater that one half, or nuclei in solid or liquid-crystal 
samples, have also been studied, but for simplicity in this 
paper I will only consider systems of spin-half nuclei in 
the liquid state. The NMR Hamiltonian for such a sys- 
tem is in general given by 



j<k 



(1) 



although in real systems the coupling strengths LUjk be- 
tween distant nuclei can frequently be taken as zero. 
Spin-spin couplings are fairly weak, and so in many cases 
(including all fully heteronuclear systems) it is possible 
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to use the weak coupling approximation 
leading to the simplified Hamiltonian 



(2) 



(3) 



jk 



where the sum is now taken over spin pairs with non- 
negligible couplings. 

NMR texts and papers typically describe spin systems 
using product operator notation 2^ 2§| , which is closely 
related to but not quite identical to conventional physics 
notation [l^. Within this language a two-spin system 
would be described by the Hamiltonian 



where 
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and the factor of h has simply been dropped by choosing 
to work in angular frequency units. 



Quantum logic gates 

It is well known that in order to perform general quan- 
tum computations, it is only necessary to implement 
single-qubit gates, which change the state of a single 
qubit, and one non-trivial two-qubit gate, for which the 
final state of at least one of the two qubits involved de- 
pends on the initial states of both qubits, so that the two- 
qubit gate encodes some sort of conditional logic [sij . 

Single-qubit gates correspond to rotating a single spin 
in its own one-spin Hilbert space. For a one qubit com- 
puter, implemented using a single nuclear spin, this can 
be achieved by applying RF fields. For simplicity if of- 
ten best to consider only resonant RF fields, so that the 
rotation has the form 



exp[— i6'(/a; cos ( 
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where 9 and cf) are the pulse nutation (fiip) and phase 
angles. Rotations about axes not in the xy plane can 
be implemented as sequences of pulses. For example, 
rotations around the z axis are easily constructed with 
composite Z-rotations js^], using identities such as 



= 90„ 9^ 90- 



(7) 



where the pulse sequence is written with time running 
from left to right, so that the leftmost pulse is the first 
pulse applied. In larger spin systems it is necessary to do 
this in a qubit-selective manner; in a fully heteronuclear 
spin system qubit selection is simple, as every spin will 
be a long way from resonance with every other spin, and 



simple pulses applied on resonance can be used, but in 
homonuclear spin systems more sophisticated approaches 
are necessary [l6|. 

Next I turn to non-trivial two-qubit gates in two-spin 
(two-qubit) systems. In NMR experiments the key two- 
qubit gate is the controlled-Z gate 



(8) 
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which is easily converted to a controUed-NOT by a pair of 
Hadamard gates. Controlled-Z is symmetric between the 
two spins and it can be easily decomposed with product 
operators [s^ as 

controUcd-Z = exp[-i7r/2 {^E - I, - + 21, S^)] (9) 

where E indicates the identity matrix. All four terms 
commute, and so can be considered individually. The 
^E term is just a global phase, and can be ignored as 
usual. Terms in I, and Sz are single qubit rotations, and 

ubit gates, or simply 
il Isil. This leaves 



so can be implemented with single-c 
absorbed into the reference frame 
the only important term, which corresponds to evolution 
under the spin-spin coupling term, ■kJ2IzSz for a time 
1/2 J. The spin Hamiltonian will include both Zeeman 
and coupling terms, but conventional spin-echo sequences 
ISOl can be used to remove the undesirable terms ISTl. 



COMPOSITE PULSES 

Composite pulses i21j have found widespread use in 
conventional NMR experiments to reduce the effects of a 
wide range of experimental imperfections, most notably 
off-resonance effects, which arise when the RF field is not 
quite resonant with the transition so that nutation occurs 
around a tilted axis, and pulse length errors arising from 
variations in the strength of the RF field. (These er- 
rors could be better described as pulse strength errors, 
but the unhelpful name is almost universal.) As similar 
imperfections are likely to affect most experimental im- 
plementations of quantum information processing there 
has been interest in applying these ideas. 

Composite pulses are only one of a whole range of tech- 
niques for improving the quality of RF excitation in NMR 
experiments. Traditionally these can be divided into 
composite pulses, made up of a small number of pulses, 
each with the same frequency and strength but differ- 
ing in length and phase, and shaped pulses [sl], which 
contain a very large number of elements, each with the 
same frequency and length but differing in strength and 
phase. A further pragmatic distinction is that composite 
pulses can often be derived and explained using analyt- 
ical approaches, while shaped pulses are frequently de- 
rived numerically using analytic ideas only as an outline 
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guide. These distinctions have increasingly broken down 
with the advent of strongly modulated composite pulses 
[3^ : these contain a small number of pulse elements, 
but these pulses are permitted to vary in frequency and 
amplitude as well as phase, and the pulse sequence is 
obtained by numerical optimisation. More recently still 
these have largely been superseded by arbitrary shaped 
pulses developed using optimal control theory, usually 
based on the GRAPE algorithm [iO, EH . In this review 
I will confine myself to "conventional" composite pulses, 
with analytical derivations. 

Most conventional composite pulses are not suitable 
for use in quantum computers, as they are optimised for 
particular classes of initial state: for example, most com- 
posite 180° pulses are optimised either for inverting the 
population of a spin or for producing a spin echo. By con- 
trast, pulses used on quantum computers must be general 
rotors, which perform well for any initial state. Com- 
posite pulses of this kind are rarely used in conventional 
NMR, but a small number of so-called Class A composite 
pulses [21I are known, and these have been developed for 
use in quantum computation. A method for construct- 
ing general rotors from conventional point-to-point pulses 
has also been described [421. 



Fidelity Measures 

The quality of a composite pulse for quantum com- 
puting can be assessed in various ways, but in practice 
there are two important families of approaches. The most 
direct approach is to expand the propagator for the com- 
posite pulse as a power series in the size of the error, and 
determine the size and order of the lowest order error 
term. As an example consider implementing a 180° ro- 
tation using a naive pulse with a fractional pulse length 
error of e, so that the flip angle of the pulse is in fact 
180 X (1 -t- e). The ideal propagator is then 



U = exp[— i7rcra;/2] 



while the actual propagator is 
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V = exp[-i7r(l + e)cra;/2 

-i i " I 7r/2 



and so the naive pulse has an error of order e. Alter- 
natively the quality can be assessed by calculating the 
propagator fldelity between U and V, given by 



F = \Tt{VU-^)\/Tt(UU-^), 



(13) 



and then expanding the fidelity as a power series in the 
error. For the naive pulse considered above the fidelity is 



and the naive pulse has infidelity of order e^. The dif- 
ference between these two methods of assessing a pulse 
must be borne in mind when comparing pulses in differ- 
ent papers; in general an error of order n will correspond 
to an infidelity of order 2n 43|, |44 1 . 



Off-resonance errors 

An early composite 90° pulse tackling off-resonance er- 
rors was described by Tycko replacing a 90° pulse 
with the three pulse sequence 385°320^^25°. This has 
subsequently been generalised to give the CORPSE family 
of composite pulses 23, [l^l, in which a 6^ pulse is re- 
placed by three pulses, applied along the +x, —x and +x 
axes as before, with flip angles given by 



V , I sm( 

9i = 2ni7r H arcsin , 

2 V 2 / 



62 = 2n27r — 2 arcsin 



sin(6i/2)\ 



S /sin(6'/2)\ 
% = 2n37r + - - arcsm ( j 



(15) 
(16) 
(17) 



F = 1 - e%78 + 0(e*) 



(14) 



where ni, n2 and are integers, with the best results 
[23 | occurring for ni = n2 = 1 and 77.3 = 0. These se- 
quences have been demonstrated by NMR [i^l, SQUID 
[46j and neutral atom 47| experiments. Pulse sequences 
have also been desigried which are tailored for particular 
off- resonance effects i48il . 



Pulse length errors 

While off-resonance errors are important in conven- 
tional NMR, they can largely be avoided in quantum in- 
formation processing experiments. Pulse length errors, 
however, remain a universal problem, arising from in- 
homogeneities in the RF field, either in space (over a 
macroscopic sample) or in time (due to slow fluctuations 
in amplifier power). There has, therefore, been consid- 
erable interest in composite pulses to tackle pulse length 
errors, which can largely be traced back to a three pulse 
composite 180° pulse due to Tycko [45| or to the BBl 
family of sequences discovered by Wimperis [i^ . Tycko's 
pulse sequence has been generalised to give the SCROFU- 
LOUS family of composite pulses but using the BBl 
family is preferable in most cases. 

BBl differs from many other composite pulses in that 
it seeks to design an error-correcting pulse, which can be 
combined with the naive error-prone pulse to give a more 
accurate compound pulse, much as a contact lens can be 
used to correct eyesight. Originally [4§| this error cor- 
recting sequence (sometimes called a Wl sequence) was 
placed before the naive pulse, but it can instead be placed 
after the naive pulse, or indeed in the middle of it [2J,[2^. 
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It comprises three pulses, in the form 180^^ 360^^ ^^^^i' 
with <^2 = 3 and 



(18) 



6i = ± arccos — — 
in 



where the choice of sign is unimportant as long as it is 
made consistently. 



In addition to NMR experiments [25|, l50| BBl pulses 
have been demonstrated in electron spin resonance [51 



in other fields 5^, 54 



55|. 



52| and neutral atoms 47 1, and have inspired applications 



Two qubit gates 

The method can be extended to build two-qubit gates 
which are robust to variations in the size of the underly- 
ing scalar coupling jsil . [60l | using an analogy between ro- 
tations on the Bloch sphere and the rotations in a multi- 
qubit Hilbert space which correspond to evolution under 
the spin-spin coupling. In combination with the single 
qubit gates described previously these provide a univer- 
sal set of robust quantum logic gates. They have been 
demonstrated using NMR techniques [i^, but it is not 
yet clear how important they will prove in practice. 



Higher precision sequences 

BBl has proved a remarkably successful composite 
pulse, and is surprisingly difficult to improve upon. BBl 
pulses can be derived by designing composite pulses 
which suppress first order pulse length errors, but it turns 
out that BBl also suppresses second order errors auto- 
matically, leaving only third order errors (sixth order in- 
fidelity). It is not clear why this fortuitous double can- 
cellation occurs, and it is not a general feature of com- 
posite puls es. Other pulses with similar properties are 
known [56], but these have no advantages over BBl. Be- 
yond this, BBl puls es are also relatively robust to off 
resonance-errors |2J|, and generally insensitive to small 
errors in their implementation, so that BBl pulses work 
in practice very much as expected from theory 25 1. 

Although BBl has proved highly successful, it is ob- 
viously interesting to seek still better pulse sequences. 



and Brown et al. have tackled this in two ways |43l . 
Firstly they have shown how the BBl approach can, in 
effect, be nested, creating ever higher orders of simulta- 
neous correction. A robust 90° pulse from the B4 fam- 
ily of pulses (which remove the third order error term) 
has been implemented in NMR experiments 25| , but this 
composite pulse is very long (the correction sequence con- 
tains 27 pulses with a total length equivalent to a 7200° 
rotation) and does not perform much better than BBl. 
Secondly they have described a general method, using in- 
sights from the Solovay-Kitaev theorem 57[ , to show how 
arbitrarily accurate composite pulses can be constructed 
in general, by building a series of correction sequences 
which correct errors one order at a time. An expanded 
version of part of their method written in more conven- 
tional NMR notation is also available 58|. Once again 



high order corrections sequences developed using these 
ideas can become extremely long, and it is not clear how 
well such complex pulses will work in practice. In com- 
parison BBl pulses are extremely robust jSH], and so may 
prove the ideal compromise between theoretical precision 
and practical implementation. 



Conclusions 

Composite pulse techniques adapted from conventional 
NMR experiments have already proved to be extremely 
useful in NMR quantum computation, and have begun 
to find wider applications in related fields. Pulse length 
errors, which ultimately relate to uncertainties in the 
strength of the external control fields, are likely to have 
analogies in many implementations of quantum informa- 
tion processing, and the BBl pulse sequence provides 
an apparently ideal method for tackling these, combin- 
ing good suppression of errors, relative simplicity, and an 
apparent robustness to imperfections in its implementa- 
tion. 

I thank the UK EPSRC for financial support. 
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